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Numerical calculations play an important role in clar- 
ifying the features in strongly correlated electron sys- 
tems. 1 In particular, quantum Monte Carlo (QMC) and 
density matrix renormalization group (DMRG) methods 
are often adopted to quantum spin and electron systems. 
However, most of them are useful only for the system 
without geometrical frustration or in low dimension. 

Exact diagonalization (ED) is applicable even in two- 
or three-dimensional frustrated systems, although there 
exist limitations to system sizes, e.g., 36 sites in two- 
dimensional spin-1/2 kagome antifcrromagnctic Heisen- 
berg systems at present. To increase the system size, cal- 
culations with some approximation have been attempted. 
One of the successful methods is a numerical contractor 
renormalization method. 2 ' 3 In this method, an effective 
Hamiltonian is constructed on a truncated local basis set 
with ED and the calculated system size increases up to 
48 sites in a spin-1/2 kagome system. 3 

In this paper, we carry out a new technique named 
projection ED and indicate that it is applicable to the 
frustrated models in the strong coupling limit with highly 
degenerate ground states in a general way. Let us con- 
sider the Hamiltonian 
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H =Ho + Ti\ ■ 



(1) 



Here, the Hamiltonian Ho has classical degrees of free- 
dom and its eigenvalues are discrete. In addition, its 
ground state is assumed to be macroscopically degen- 
erate and to have an excitation gap A (see Fig. 1(a)). 
Hi is a quantum term that hybridizes the eigenstates of 
Ho, which will be treated perturbatively. 

As typical examples, we consider frustrated XXZ spin- 
1/2 Heisenberg and spinless fcrmion t-V models in a 



strong coupling region, 



\Jxv\ ^ Jz m the former 



model and \t\ <C V in the latter model, where J xy and 
t give perturbation terms Hi. The term Ho corresponds 
to an Ising model. On some geometrically frustrated lat- 
tices, i.e., on triangular, kagome and pyrochlore lattices, 
the models with Hi = have macroscopically degenerate 
ground states when interactions are antiferromagnetic, 
i.e., J z > or V > 0. 4 ~ 7 

For A^> Ei, where Ei is the energy scale of Hi , we may 
apply the perturbation theory to describe low energy 
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Fig. 1. (a) Macroscopic degenerate ground state with excitation 
energy A. (b) Calculated energy within the projected subspacc 
for ei < A. 



states as those spanned within the Hilbcrt space pro- 
jected to the degenerate ground states of Ho, as shown 
in Fig. 1(b). Within the second order of the perturbation 
of Hi, we can obtain an effective Hamiltonian: 

H=E + VH 1 V -VHxQ—QHiP, (2) 

Ho 

where V is a projection operator, which takes 1 in case 
that the applied state is one of the ground states of the 
Ising model, and for higher energy states. Q is the 
complement of V and satisfies V + Q = 1. Eq is the 
ground state energy of the classical model Ho- 

In these cases, we have difficulties caused by (i) macro- 
scopic degeneracy in the ground state and (ii) the nonlo- 
cal property of V . Due to the former reason, analytical 
approaches for perturbative calculations are extremely 
difficult. Therefore, we need to resort to numerical cal- 
culations. To the advantage of ED, the number of ground 
states of Ho is much smaller than the total dimension of 
the original Hilbcrt space of H, in general. Thus, it be- 
comes possible to handle a much larger system than in 
the case of conventional ED. 

The latter difficulty originates from the fact that, un- 
like the simple cases such as the Hubbard model at 
U — ► oo, the unperturbed ground states themselves are 
highly nontrivial. This gives complex forms for analyti- 
cal expressions of the projection operators V and Q, and 
hence, the overall perturbation terms. In addition, since 
the ground states depends on lattice structures, the for- 
mula also depends on them and is thus nonuniversal. 

To avoid such difficulties, we introduce ED, which con- 
sequently makes it possible to treat the problem in gen- 
eral. Hereafter, we show the details based on the calcula- 
tions for the spinless fermion t-V model on a triangular 
lattice as an example. A Hamiltonian can be written as 



Ho =V 



(ij) 



c\Cj +H.c 



(3) 



(4) 



Here, c\ (cj) is the creation (annihilation) operator of 
fermions at site i and n, = c\ci. V is a repulsive interac- 
tion between nearest neighboring sites. We allow hopping 
between nearest neighboring sites on a triangular lattice 
with an amplitude t > 0. 

Lanczos procedure has been performed by applying the 
Hamiltonian Hi, which is an O(M) calculation. Here, M 
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N 


Al (ED) 


M (projection ED) 


12 


924 


66 


18 


48,620 


386 


24 


2,704,156 


2,142 


30 


155,117,520 


13,694 


36 


9,075,135,300 


84,486 


42 


538,257,874,440 


553,574 


48 


32,247,603,683,100 


3,720,042 



Table I. Number of states at half-filling for conventional and pro- 
jection EDs on N-site clusters for the triangular t-V model. 
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Fig. 2. (Color online) Energies E/tN at density p on triangular 
spinless fermion t-V model for V/t = 50. The results for up to 
48-site clusters are shown. 

is the number of ground states of TLq at a given fermion 
density. For the hopping from site i to j, a locally mod- 
ified energy can be obtained by confirming the 10-site 
configuration including 8 neighboring sites of i and/or j 
on a triangular lattice. When the modified state belongs 
to the ground state of TLq, the hopping term contributes 
as the first-order perturbation term in eq. (2). In the case 
that the modified state is an excited state of Ho , further 
hopping processes which return the excited state to one 
of the ground states, lead to the second-order perturba- 
tion term in eq. (2). Generally, the possibility of such a 
correlated hopping can be confirmed by considering the 
13-sitc configuration at most. Since what we need is the 
energy of a locally modified state for the Ising model, the 
type of calculation is an 0(1) calculation. In this manner, 
projection ED is generally an O(M) calculation. Note 
that, even if we got the formula in eq. (2) analytically, 
its direct Lanczos operation is also O(M). 

As shown in Table I, the number of states M is much 
smaller than that of conventional ED. Thus, we can reach 
a much larger cluster than in the case of conventional 
ED. For example, we could carry out the calculation up 
to 48-site clusters in the second order of perturbation 
easily. Energies E/t as a function of fermion densities 
are shown in Fig. 2. Here, the calculations have been 
performed for the numbers of sites N = 18, 24, 30, 36, 42 
and 48 at the densities 1/3 < p < 2/3 for V/t = 50. Note 
that we can handle a much larger system on kagomc and 
pyrochlore lattices by also applying projection ED. 

The accuracy of our method can be confirmed by com- 
parison with conventional ED in small clusters. We calcu- 
late the energy differences between projection and con- 
ventional EDs using 8 = \(E — E)/ E\, where E is the 
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Fig. 3. Comparison between conventional and projection EDs on 
the spinless fermion t-V and hard core boson models on the 
triangular lattice. The calculation has been performed at half- 
filling for both first- and second-order perturbations. For V/t > 
10, the energy differences are proportional to (t/V) 2 for the first- 
order calculation and (t/V) 3 for the second-order calculation. 
The fitting results are shown by the lines. 

energy calculated by projection ED, and E is the en- 
ergy obtained by ED. The results for a 24-site cluster 
are shown in Fig. 3. In the parameter region, where the 
n-th order of perturbation is sufficiently accurate, E is 
consistent with E within 0((t/V) n ), i.e., the errors 8 are 
proportional to (t/V) n+1 . As shown in the figure, we can 
obtain a good agreement in the second-order perturba- 
tion for V/t > 10, where the errors are proportional to 
(t/V) 3 . In the figure, the results for a hard core boson 
model on a triangular lattice are also shown. 

Finally, we mention that researches of the hard core 
boson and t-V models on the triangular lattice by pro- 
jection ED are significant. These models have been paid 
much attention to, after the prediction of the realization 
of a new type of ground state, i.e., a supersolid state in 
the hard core boson model 8-10 and a pin-ball liquid state 
in the t-V model. 11 In the former, most of the calcula- 
tions have been performed by QMC simulation, and in 
the latter, ED has been performed. In QMC, it is diffi- 
cult to calculate in parameter ranges V/t > 10. For ED, 
the system size is limited. Thus, our proposal should be 
important to clarifying the features at certain parame- 
ter ranges or system sizes, at which previous methods 
cannot be applied. 
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